charpits method pde